Shot noise in diffusive conductors: A quantitative analysis of electron-phonon 

interaction effects 
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Using the 'drift-diffusion-Langevin' equation, we have quantitatively analyzed the effects of elec- 
tron energy relaxation via their interaction with phonons, generally in presence of electron-electron 
interaction, on shot noise in diffusive conductors. We have found that the noise power Si(ui) (both 
at low and high observation frequencies u>) drops to half of its 'mesoscopic' value only at (5 > fOO, 
where (3 is the ratio of the sample length L to the energy relaxation length Z p h (the latter may be 
much larger then the dephasing length), ft means in particular that at low temperatures the shot 
noise may be substantial even when L ~ fO" 2 - fO" 1 cm, and the conductor is 'macroscopic' in any 
other respect. 
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come important at L > l ee = y/Dr t 
intensity only slightly at low frequencie; 



Macroscopic metallic conductors, in which electrons 
are in local equilibrium, do not exhibit non-equilibrium 
("shot") noise at low voltages. On the other hand, short 
diffusive conductors (with length L much shorter than 
the effective electron-phonon relaxation length l p h — 
y/DTph) do show this effect as soon as the average en- 
ergy eV/2 acquired from the electric field becomes larger 
than the ener.gy scale T of equilibrium thermodynamic 
fluctuationscTEl. Electron-electron scattering (which be- 

) affect the shot noise 
1, though at fre- 
quencies above the reciprocal Tiaouless time TjT 1 = D / L 2 
their effect is more substantial!! This low sensitivity of 
shot noise to electron-electron scattering, as well as to 
screeninga, may be readily understood: electron-electron 
interaction cannot drain the energy supplied by the ex- 
ternal electric field from the electron subsystem, which is 
thereby "overheated". Shot noise may be considered as 
a direct result of this deviation from equilibrium. 

The effect of electron-phonon interaction is quite dif- 
ferent: it may drain extra energy from the electron sub- 
system, bring it closer to local thermal equilibrium, and 
hence suppress the shot noise. The analysis of these 
effects has been- addressed in many recent works, both 
analyticallw jQullJ and using numerical Monte-Carlo 
simulations^ (see also Ref.lia). While describing the re- 
duction of shot noise with growing values of j3 = L/Z p h, 
these works cannot provide realistic values for the noise 
at P > 1. For example, assuming a simplifying form of 
the electron-phonon interaction, it was concluded in Ref.El 
that the low-frequency noise behaves as Sj ~ e~~@ . On 
the otherdiaxid, a crude partition of the sample into L/lp^ 
segmentsLffl'Ell, or, alternatively, the assumption of a con- 
stant relaxation timet!, lead to the conclusion that the 
noise should vanish much slower, as since such ap- 
proaches neglect the dependence of the relaxation rate 
on the excitation enertfK. This dependence was taken 
into account by Nagaev&Q. By asymptotically solving the 
Boltzmann equation, he showed that for the case of zero 
temperature and zero frequency, and in the limit (3 1, 
the noise intensity behaves as Si ~ /3~ 2 / 5 . However, the 



full dependence of the noise on /3, as well as the more real- 
istic case of finite temperatures and frequencies, has not 
yet been studied. Such a study was the goal of this wocJ< 
We have used the 'drift-diffusion-Langevin' equation! 
to calculate the current noise power as a function of the 
ratio (3, the temperature T, and u>, for the two limiting 
values of the ratio 7 = L/l ee . __ 

The drift-diffusion-Langevin equatioroB is based on 
integration of the Boltzmann-Langevin equationlij over 
the electron energies. Its self-consistent solution with 
the Poisson equation that accounts for screening in the 
system enables one to calculate noise power at arbitrary 
frequenciesp. The outcome of this approach may be sum- 
marized as the following simple recipe: for a conductor 
with uniform cross-section, the spectral density of current 
fluctuations at any point x (along conductor's length) is 

OA rh 

Sr{x; u>, T) = / \K(x, x'-uo)\ 2 S{x'- T) dx' (1) 

with A the cross-sectional area. The local noise correla- 
tor S{x; T) is given by 
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a(x) dEf s (E,x;T)[l-f s (E,x;T)} (2) 
Jo 



S(x;T) = 2 



with a(x) the local conductivity and f s (E,x;T) the 
momentum-symmetric part of the local (steady state) 
distribution function at a total energy E. The response 
function K(x,x';u>) equals 1 at zero frequency, but at 
finite frequencies it is dependent on the specific geome- 
try of the conductor and its electrodynamic environment, 
always obeying the following sum rule: 



L 



K(x, x '; u>) dx' = 1. 



(3) 



In this work, we concentrate on the geometry which 
looks most promising for experimental observation of the 
considered effects, namely, a thin and long conductor 
close to a ground plane. For the quantitative applica- 
bility of our results, thickness of the conductor and its 
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distance from the ground plane should be much smaller 
than L, but its second transversal dimension (parallel to 
the ground plane) may be arbitrary. (As a result the 
model is applicable, e.g., to a two dimensional electron 
gas gated by a close electrode, which also serves as the 
ground plane) . In this geometry the response function £pr 
noise current in the external electrodes is very simpleQ: 



K e {x'\uj) = K 



L cosh(K2;') 
2"sinh(Ki/2)' 



(4) 



Here re = yJ—iu/D', D' = D + aA/C , C is the (dimen- 
sionless) linear capacitance between the conductor and 
the ground plane, and D is the diffusion coefficient. The 
non-equilibrium distribution function f s (E,x) should be 
found by solving the stationary Boltzmann equation. At 
this stage it is convenient to use dimensionless quantities 
£ = x/L, e = E/eV, and t = T/eV. In the diffusion 
limit (i « L), the Boltzmann equation is 



(D/L 



d£ 2 



(5) 



with I(e, £) the collision integral. 

We study here two limiting cases. In the first limit, 
electron-electron scattering within the conductor is neg- 
ligible, 7 = L/l ee -c I , so the collision integral involves 
only scattering by phonons. For the deformation poten- 
tial scattering it may be presented aaifl 

J( £ ,0 = — / duu? {(1 - /,) [/+(1 + N) + f-N] 
T v Jo 

-f s [(l-f-)(l + N) + (l-f+)N]} (6) 
/.(eiw.O, N = N(lo) = 



with U = /.(e,0, ft = 
l/[exp(u/t) — I], and 
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27T \2hkpv s 



Ti 2 pv s 



(7) 



with kp the Fermi wavenumber, m the effective mass, S 
the deformation potential constant, v s the sound velocity, 
and p the mass density of the material. Since the voltage 
drops only across the conductor, the distribution func- 
tion must approach the equilibrium distribution at the 
conductor-electrode interfaces. The boundary conditions 
for Eq. (§) are therefore 



/ s (e,Tl/2) = /o(£Tl/2) 
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1 + exp 



(8) 



In the second limit, electron-electron scattering is 
strong, 7 ^> 1, so the electrons are locall*i-|thermalized. 
The distribution function is then given bytJ 



1 



1 + exp 



t e (x) 



(9) 



The equation for the electron temperature t e (x) = 
T e {x)/eV in this limit was obtained in Refo by multi- 
plying the Boltzmann equation by e and integrating it 
with respect to energy. The resulting equation reads 
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where C(5) — 1.04 is the Riemann zeta function. 

Ty which appears in equations 10) is the energy re- 
laxation time of electrons with typical energy eV. The 
relaxation length is then / p h = y 1 Dtv ■ One can imme- 
diately notice from equations (||-§|) an d (p[-p^|) that the 
dependence of f s (E,x) on the physical variables of the 
problem eV, T, L, and l p h in each of our cases is oi 
through the parameters t = T/eV and (3 — L/l p ^. 
From equations it is then seen that for a uni- 

form conductor [cr(x) = a] the only additional parameter 
which affects the normalized noise value a = Si/2el is 
\k\L = ^Jujtti with tt = L 2 /D' the effective Thouless 
time. In particular, at fixed ujtt the only dependence of 
the noise on the sample length is due to its ratio (3 with 
the thermalization length. 

Equation (||) can be solved analytically in three limit- 
ing cases: (3 — 7 = (no energy relaxation), (3 — 0, 7 — > 
00 (strong local thermalization), and (3 — > 00 (strong 
thermalization to lattice temperature). The correspond- 
ing distribution functions are 

/.(e,0 = (l/2 + 0/o(e+l/2) 

+ (l/2-0/o(e-l/2) /? = 7 = 0, (Ha) 



1 + exp 



1 + exp 



t h (0 



(3 = 0, 7^00, (lib) 



t 



(3 -> 00, 



(11c) 



with 

Vi 2 + 3(l-4£ 2 )/4vr 2 



the hot-electron temperature th(£) = 
The frequency and temperature 



dependences rof the noise power in these limiting cases 
were given iria. 

In order to study the crossover region {i.e., at finite 
values of 0) we solve equations (^|) and ( |l0| ) numerically. 
Results for the distribution functions for several values of 
(3 are presented in figures 1 and 2. Figure 1 shows the de- 
pendence of f s {^,0 on energy and position at t = 0.01, 
with no electron-electron scattering. Plots (a) an d (f ) 
corre spon d to the limiting cases described by Eqs. ( |lla| ) 
and (lie), respectively. At large (3 the integral in Eq. (|6J) 
must be small, and so the distribution is close to a Fermi 
Dirac distribution in most of the sample. However, the 
effective temperature of this distribution approaches the 
lattice temperature t = 0.01 only at (3 > 10 5 . This form 
of the distribution function means that a perturbation 
expansion of Eq. ([5]) around the distribution at {3 — * 00 
gives results which are far from reality for values of (3 
smaller than ~ 10 5 . Note also that at distances of the 



2 



order L//3 1 / 2 from the edges of the sample the distribu- 
tion is much sharper than in the bulk. We will show that 
this fact has strong implications on the high frequency 
noise. 



y«1 






Figure 2(a) shows the distribution function for the case 
of strong electron-electron scattering, and (3 = 10. At any 
other (3, the plot looks generally the same, with the only 
difference being that the width of the distribution on the 
diagonal e + £ = would change according to the local 
electron temperature t e (£;(3). The dependence of this 
temperature on position and on (3 is shown in Fig 2(b). 

Figure 3 shows our main result: the normalized spec- 
tral density as a function of (3 for various frequencies. 
ujTT,t <C 1 and (3 ^> 1, Nagaev's asymptotic results' 
are reproduced; in our notation they read 



Si/2el : 
Si/2el - 
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(7 « 1), 
(7 » 1). 



(12a) 
(12b) 
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However, for the case 7 < 1 the result (12a) relies on 
the fact that at large (3 the electron distribution function 
is invariant to diagonal transformations £ — > £ + £0 , £ 
e — £o- As mentioned above [cf. Fig. 1(e)], this invari- 
ance is not valid at distances £ ~ /3 _1//2 from the edges of 
the sample. Due to the form of the response function (Q) 
the high frequency noise is sensitive to the distribution of 
electrons only at distances £ ~ 1/\kL\ — (ujtt) -1 ^ 2 from 
the edges. Therefore, the noise behaves with its asymp- 
totic power-law form only at f3 3> wtt, as clearly seen in 
Fig. 3 (solid lines). 



FIG. 1. The symmetric part of the electron distribution 
function f s (s, f) for different values of /3 = L/l p h and for 
weak electron-electron interaction (7 = L/l ee <C 1). For all 
curves, T = O.OleV. 
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FIG. 2. (a) Same as in Fig. 1 but for the case 7 >• 1. (b) 
The dependence of the electron temperature t e on the position 
and on the electron-phonon interaction parameter f3. 
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FIG. 3. The dependence of the noise spectral density on /3 
for t = T/eV — 0.01 at various frequencies, for the two cases: 
L <C lee (solid lines) and L 3> lee (dashed lines). The same 
dependence for t — 0.03 and zero frequency is also shown. 
Inset: the dependence of the spectral density on frequency at 
t = 0.01 and = 100. 

A similar analysis can be done for the case 7 ^> 1 
(dashed lines in Fig. 3). Here Si ~ /? -2 / 5 because / s (e, £) 
does not depend on £ at j3 ^S> 1 . As can be deduced from 
Eq. (|l0|), this constancy of f s does not hold at distances 
smaller than £ ~ /?~ 2 / 5 from the boundaries, and so, by 
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the same arguments a s for 7 < 1, the high frequency 
noise behaves as Eq. (|l2b|) only at /? > (wt t ) 5/4 . In 
contrast to the case of 7 1, however, in this case the 
noise at 1 <C (3 < (cjt<t) 5 / 4 does not retain its /3 = 
value, because the electron temperature decreases with 
(3 also near the boundary, see Fig. 2(b). The inset of 
Fig. 3 shows the frequency-dependence of the noise spec- 
tral density at /3 = 100 and t = 0.01. More details of the 
dependence of the noise on frequency and temperature 
at P = were given elsewherefl. 

To summarize, we have performed numerical calcu- 
lations of the non-equilibrium noise in samples with 
strong elastic scattering and with an arbitrary strength 
of electron-acoustic phonon scattering. We have shown 
that the shot noise decreases very slowly with the sample 
length, and is of the order of the Schottky value even for 
L/lph ~ 100. In experiment, it means that when deal- 
ing with otherwise 'macroscopic' samples (at low tem- 
peratures, of the order of one millimeter), the a priori 
assumption of vanishing shot noise may be wrong. In 
view of the persistent improvements in accuracy of noise 
measurements in recent yearal3~Eil and the growing gen- 
eral interest in the relationSj-hetween inelastic relaxation, 
fluctuations, and dephasingE3, we believe that there is a 
considerable chance of experimental confirmation of this 
prediction. 

The work was supported in part by DOE's Grant 
#DE-FG02-95ER14575. 
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